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ABSTRACT 

The analysis presented in this paper verifies which of the mixing rules are 
best for real components of interstellar dust in possible wide range of wave- 
lengths. The DDA method with elements of different components with various 
volume fractions has been used. We have considered 6 materials: ice, amor- 
phous carbon, graphite, SiC, silicates and iron, and the following mixing rules: 
Maxwell-Garnett, Bruggeman, Looyenga, Hanay and Lichtenecker which must 
satisfy rigorous bounds. The porous materials have also been considered. We 
have assumed si mplified spatia l distribution, shape and size of inclusions. The 



criteria given bv lDraind (|l988f ) have been used to determine the range of wave- 
lengths for the considered mixtures in order to calculate the Qoxt using the 
DDA. From all chosen mixing rules for the examined materials in majority 
of cases (13 out of 20) the best results have been obtained using the Lichte- 
necker mixing rule. In 5 cases this rule is better for some volume fraction of 
inclusions. 

Key words: ISM: general, dust, extinction, interstellar grains, mixing rules, 
discrete dipole approximation 



1 INTRODUCTION 

Interstellar dust is a mixture of grains of different shape, size and chemical composition. 
The most frequent are assumed to be different allotropic forms of carbon, a-SiC, "astro- 
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nomical silicate" and ice. Many authors proposed grains which have been mixtures of vari- 
ous materials. In order to obtain the refractive indices of such types of grains the formulae 
describing different mixing rules are used. Those formulae have been derived for various 
assumptions concerning arrangement of inclusions and their shape. The best known exam- 
ples of the effective medium theories (EMT) are theories by Maxwell- Garnett and Brugge- 
man. The Maxwell-Garnett mixin g rule has been re-derived by Bohren and Wickramasinghe 
(jBohren fc Wickramasinghd 119771 ) for the spherical inclusions arranged chaotically. The as- 
sumption of spherical shape of inclusions or generally separated inclusion structure causes 
asymmetry of this mixing rule, whereas the formula derived by Bruggeman (and a similar one 
by Landauer) with non-spher ical inclusions, tightly adjo ining to each other aggregate struc- 
tures leads to full symmetry (IChylek &: Srivastavalll983r). Also in deriy i ng the Looy e nga and 



Hanay rules particular models of mixtures were used (ILoovenga 



19651), 



Beek f 



19671^ 



Sihvola 



19731 )). Many existing mixing rules were described by iBeekl (119671 ) or ISihvolal (1l973l ). for 
example. The Lichtenecker mixing rule has one drawback that it was derived on the basis of 
fitting to the empirical data. It lacks a physical model apart from some theoreti cal justifica- 



tion b ound with an artificial decomposition of geometrical shapes of inclusions (jZakri et al. 



(119981 )). It is worth mentioning that the interaction between inclusions for small volume 
fractions is limited and may be omitted. On the other hand for large volume fractions of 
inclusions this effect is not negligible. Those interactions or their lack are included in more or 
less explicit way in the assumptions for the mixing rules. Therefore, they lead to the limited 
applicability of a particular mix i ng ru le. The influence of interactions between inclusions 



was discussed by iPerrin &: Lamyl (Il990l ) for Maxw ell- G arnet t and B ruggeman mixing rules. 
In case of derivation of Looyenga rule the author (ILooyengal (119651 )) avoided the discussion 
on those interactions. Most of the mixing rules may be applied with good approximation 
for various components with small volume fraction of inclusions (up to a few per cent) 
when the average distances between inclusions are large and the interactions are small. For 
higher volume fractions the choice of a mixing rule is difficult and very important. We used 
the DDA method with elements (dipoles) of different components (refractive indices) with 



various volume fractions in order to ca^ 



authors (IChvlek et all mm 



lati et al 



culate t he extinction coefficients fo r grains. Many 



(I2OOJ) , IVoshchinnikov et al.l (120071 ) ) compared the 



extinction coefficients calculated in this way with the extinction coefficients calculated from 
Mie theory for grains of refractive indices computed using various mixing rules in order to 
choose the best rule. The extinction coefficients obtained from the Mie theory are the same 
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as those calculated with the DDA method using a large number of dipoles. Using a large 
number of dipoles decreases the influence of granularity of grains but at the same time it 
requires a very long computing time. In order to avoid the influence of the method for cal- 
culating the extinction and shorten the c omputing tim e we used the DDA method for both 
cases. Many authors (cf. Beekl ( 1967 ) and Sihvola ( 1973 )) compared experimentally obtained 
refractive indices of mixtures with results given by various mixing rules. Depending on the 
volume fractions of inclusions, their kinds, shape, spatial distribution and frequency range 
different mixing rules fitted experimental permittivity. Our choice of mixing rules is based on 
numerical experiment which allowed us to examine the mixtures in a wide frequency range. 
In this paper we have assumed random spatial distribution, pseudospherical shape and one 
size of inclusions. However, in this numerical experiment usi ng the DDA meth od and the 
Rayleigh and non-Rayleigh cluster dipol inclusions method (jWolff et al.l (119941 )) we could 
change the geometrical parameters of inclusions and their spatial distribution according to 
assumptions leading to different mixing rules and verify their applicability. 



2 PREPARATION OF OPTICAL DATA 

In this work 6 materials have been considered: ice, amorphous carbon, g raphi te, SiC, silicates 



and iron. The refractive indices of ice have been taken from 
carbon from 



Zubko et al 



WarrenI (119841 ). of amorphous 



(119961 ) . for graphite, SiC and "astronomical silicate" from Draine 
(|http : / / www. astro . princeton . edu / ^dr aine / dust / dust . diel. html | ). The refractive indices of 
iron have been compiled by lLynch fc Hunted (1199 ll ). We have verified whether the Kramers- 
Kronig relation is fulfilled. The differences of values for n taken from the cited literature, 
except for iron, and those calculated from Kramers-Kronig relation are negligibly small. For 
iron there are significant differences due to different methods used by various authors for 
different wavelength ranges. In order to obtain a homogeneou s data sets of refractive indices 



the values of n have been taken from 



Lynch &: Hunterl (Il99ll ) and used to calculate k values 



from Kramers-Kronig relation. We have interpolated 100 values of n and k in the range from 
0.0443 -7- 150/im and therefore obtained the same wavelength range for further calculations. 



3 MIXING RULES FOR TWO CONSTITUENTS 

Our study has been limited to grains without electric charge, magnetic suscept ibility and only 



two component mixtures. We have studied 5 mixing rules described in detail in 
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(I2OO5I ) excluding the Rayleigh mixing rule mo dified by iMeredith fc Tobiaa (Il960l ) because 
this rule has been derived for ordered mixtures (jSihvolalll973l ). The following rules have been 
taken into account: 
Asymmetrical 



Maxwell-Garnett (IBohren fc Huffman! Il983l ) 



Si -\- '2,Sy,i f (^j E 

Hanai-Bruggeman (called Hanai in this paper) (lBeeldll967l ) 
Symmetrical 



(1) 



(2) 



Bruggeman (IBohren fc Huffman! 1 19831 ) 



Ei + 2e 



Looyenga (ILooyenga 



0, 



19651) 



es=fef + {l-f)El, 



(3) 



(4) 



Lichtenecker (ILichteneckerl Il926l ) 



^OgE = f\ogEi+ {1- f)\ogEm. (5) 

In all formulae / is the volume fraction of inclusions and Em, Ei and e (without a sub- 
script) are the complex dielectric permittivities of a matrix, inclusion and mixture, respec- 
tively. 



4 WIENER AND HASHIN-SHTRIKMAN BOUNDS FOR THE 
EFFECTIVE COMPLEX PERMITTIVITIES 

Mixing rules must satisfy rigorous bounds which for complex permit ti vities o f comp osite 
of two isotrop ic components have been generalised by BergmarJ ( jl98ol ). iMilton (1980) and 
Aspend ( 1l982l ). It is necessary to discuss three cases of bounds: 

1. If we do not know the volume fraction of components and the micro-structure then 
the resulting permittivity of a mixture is located on a complex surface limited by Wiener 
bounds: 

(a) when there is no screening (all borders of inclusions are parallel to the external 
electric field) 
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matrix 



Figure 1. The area where the resulting permittivity is located for case 1. 



£p = f£a + (1 - /) Sb 
£p = fSi + (1 - /) Bm 



(6) 
(7) 



(b) with maximum screening (all inclusion borders are perpendicular to the external 
electric field): 

(1-/) 



1 = / 



1 _ / I (1-/) 



(8) 
(9) 



2. If we know the volume fraction / of the mixture components and their permittivities 
and eb the resulting complex permittivity is located in a smaller area Q'. The area Q' is defined 
by arcs of circles crossing the three points of ep{f), es(/) and ea or eb- The permittivities of 
the mixture components and e?, may have the following values: 



— ^i, — Cm or ea 



where and are permittivities of inclusion and matrix, respectively. The above rigorous 
bounds are called Hashin-Shtrikman bounds. 

3. If the micro-structure is known it is possible to further limit the area in which the 
resulting permittivity of a mixture must be located. 

In this paper the case (2.) has been considered because the volume fraction of inclusions 
is known but the information about the micro-structure is not available for some mixing 
rules. 
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matrix 



Figure 2. The area where the resulting permittivity is located for case 2. 



5 CRITERIA OF DDA APPLICATION 

The criteria for application of DDA have been described in details by 
general the criteria are as follows: 

(i) The influence of surface granularity. 
The number of dipoles N must satisfy that 

iV>iV^.nl~60|m-lp(^)-^ 

where A = 0.1 is the fractional error. 

(ii) Skin depth. 

> Nmin2 



Drainel fllQSSh . In 



(10) 



3 ' ' ^o.r 



(iii) The influence of magnetic dipole effects. 



N> N, 



min{magn) 



(90 ■ A)3/2 



eqj 



V90- A 



(11) 



(12) 



where k = ^ and acq = O.lSyum. The above equation combined with the criterion of influence 



of skin depth gives: 



N> N, 



An 



min2 



— {ka,,f\mW—)-\l + 
3 ^ ^1__SQAL 



Imp A 



(13) 



- n.1 ' ' 367r 0.1- 

The criteria given by lDraind (119881 ) (Equations [TO] and [T3i) allow to determine the conditions 
which must be satisfied in order to use the DDA method. One may either calculate the 
smallest number of dipoles at a given wavelength range or determine the wavelength range 
for a given number of dipoles. In our case the criteria have been used to determine the range 
of wavelengths for the considered mixtures in order to calculate the Qcxt using the DDA. 
For this purpose two mixture cases described by equations (6) and (8) have been used. 
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N''^^(X) eq. 10 - no screening 

^''miii2^^' cq. 13 - no screening 

^^mini^^^ eq. 1 - max. screening 

(\) eq. 1 3 - max. screening 




0.01 



0.1 J^mhl 1 '-m.il 

>-mii.2 wavelength [urn] 



Figure 3. An illustrative example showing the choice of the wavelength range which is in accordance with Draine criteria. See 
text. 

For the calculated values of ep and eg for the given volume fraction of inclusions / in 
the wavelength range from 0.0443 to 150/im the minimum number of dipoles A^mmi(-^) 
^miniW have been calculated from (fTOj) and A^mm2(-^) ^min2W from (fT3|) . Figure [3] 
shows an example of the relations of those values for 30% inlusions of graphite in amorphous 
carbon matrix. In the wavelength range for which the following ralations are simultaneously 
satisfied the Draine criteria are also satified: 

. iV^„,i(A) < 1791 
. iV^„i(A) < 1791 
. iV^,„2(A) < 1791 
• A^™„2(A) < 1791 

The two chosen values have been the limits of DDA applicability for the given mixture in 
terms of composition and volume fraction of inclusions /. The choice of wavelength ranges 
have been applied to f = 0.05 — 0.50 with step of 0.05. The value of / has been limited to 
0.5 because it is the approximate value of percolation threshold. The same procedure was 
carried out for all mixtures although after reaching the value / = 0.5 the whole procedure 
has been carried with the roles of matrix and inclusions interchanged. 



6 DETAILS OF CALCULATIONS 



The details of calculation have been given in the previous paper ( iMaron fc Marorul2005l ) with 
the only difference that we have chosen the mixtures of amorphous carbon, graphite, SiC, 
"astronomical silicate" and iron in ice and mixtures of ice in those materials for examination. 
Besides, we have considered each material with pores containing vacuum. The calculations 
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have been carried out for volume fractions from 5% to 50% with step of 5%. Draine has 
pubhshed a new version of DDSCAT program but there were no changes in the domain that 
was interesting for this work. Therefore in current work as in the previous we have used 
the version DDSCAT 5al0. The values of permittivity for mixtures calculated according 
to equations (l)-(5) have been chosen with respect to the Wiener and Hashin-Shtrikman 
bounds (it was important for Hanay and Bruggeman rules where there were 3 or 2 solutions 
of mixing equations, respectively, and also in case of porous ice). The wavelength ranges for 
each mixture have been bound according t o Draine cr i teria which is seen in Figures [9ll28l 



Many authors (including iDraind (119881 ) and 



Wolff et al. 



I994J )) have examined the influence 



of the number of dipoles on physical convergence of extinction coefficient. Numerical tests 
indicate that when the number of dipoles in the DDA method aproach infinity the extinction 
coefficient for the pseudosphere Q^^^ aproaches the extinction coefficient obtained from the 



Mie theory (Qext ~^ ^) = Qexf)- Of course, using a large number of dipoles in the 
grain implies inclu sion consisting of many dipoles a s well because its ra dius should not be 



Wolff et al 



(1l998l ) examining porous 



smaller than 50 A ( Bohren &: Wickramasinghd 119771 ). 
grains concluded that extinction obtained by DDA method for single dipole vacuum inclu- 
sion was in good accordance with that obtained from Mie theory and effective extinction 
coefficient calculated from effective medium theory, while for multi dipole vacuum inclu- 
sions they were strik i ngly n ot compliant with each other. Similar conclusions were made by 



Voshchinnikov et al. 



(120071 ) who stated that if the inclusions were not simple dipoles in the 
DDA terms the scattering charcteristics of aggregates were not well reproduced by the EMT 
calculations. Therefore, the inclusions have been assumed to be simple dipoles. We have 
limited the number of dipoles in the grain to a rather small number (1791) for the following 
reasons: 



(i) The size of single dipole inclusions (r = lOOA) are sufficient in order to safely use the 
bulk dielectric function. 



(ii) Tests by (lDrainelll988l ) showed that if the criteria (eq. [TO]) and (eq. [13]) are satisfied 
we obtain a good agreement between extinction coefficients from DDA with those calculated 
from Mie theory. 

(iii) In our approach we do not pursue the convergence of DDA results with those obtained 
from Mie/EMT but only the single dipole inclusions DDA with DDA/EMT. Therefore, the 
limited number of dipoles is sufficient providing that the Draine criteria are satisfied. 
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Figure 4. Schematic representation of grain (left - white circles depict matrix dipoles and the black ones dipole inclusions, 
right - grey circles are uniform dipoles) 



dipoles. 



Wolff et al. 



Wolff et al. 



(I1994I ) the inclusions indeed do not have to be 



As i t has been s tated by 

( I994I ) were the first to consider nonspherical and nondipole inclusions 



containing a large number of dipoles. This aproach is justified in case of fitting the theoretical 
extinction curve obtained from the DDA method with the observed one. Nevertheless, in 
this paper we compare the extiction curve calculated for grains composed of two kinds of 
dipoles with the one calculated for grains obtained by using a mixing rule. This situation 
is illustrated by the figure |H In Fig. [Deleft) the white circles stand for matrix dipoles (em) 
and the black ones for dipole inclusions (ei) of the mixture with 20% of inclusions. In Fig. 
m^right) the grey circles stand for uniform dipoles with permittivity obtained from mixing 
rules (e) for the same volume fraction of inclusions. Both cases have been computed using 
the DDA. Therefore, it is justified to use the idealised grains and inclusions. 

The arrangement of inclusions (DDA elements) is random. In order to obtain random 
number of DDA elements in the discrete dipole arrays we have used random number gen- 
erator " Research Randomizer" available at \http:/ / www, randomizer, orc^ We have generated 
10 series of numbers corresponding to the given volume fractions of inclusions out of all 
1791 dipoles. The generated numbers have been sorted in ascending order. The location of 
dipoles in the array for the spherical particle obtained from the routine calltarget.f was the 
same as for homogeneous grains with the only difference that for the randomly generated 
numbers of the DDA elements they had the refractive index of inclusions and the remaining 
elements had the refractive index of a matrix. Certainly the influence of inclusions topology 
on extinction might exist but the authors of the considered mixing rules have assumed a 
statistical distribution of inclusions. Therefore, in our calculations the random distributions 
have been used. Because there is a scattering of results for different random distributions 
the efficiency factors for extinction have been calculated for 10 different distributions and 
then averaged. For grains with radii of r = 0.15/im 10 values of efficiency factors for extinc- 
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tion Clfyf" depending on random location have been calculated. The calculations using the 
computer program DDSCAT.SalO have been carried out for wavelengths in the range per- 
mitted by the Draine criteria for the given refractive indices. The subscript i in the symbol 
Glfy^ corresponds to the number of random location of inclusions, j - the volume fraction of 
inclusion and the subscript 1 corresponds to wavelength. Next the mean extinction has been 

calculated as: 
1 1° 

Q\T' = YQT.Q7:-t- (14) 



We have calculated the standard deviation of the mean: 



i=iWi,j ) , 

where 1=10 is the number of random locations. 

The extinction for homogeneous grains has been calculated using DDA assuming that all 
elements (dipoles) consist of the same mixture with averaged refractive index calculated from 
the given mixing rule. In the obtained extinction coefficient for the homogeneous grain Qf°™°^ 
the subscripts 1 and j denote the wavelength and volume fraction of inclusion respectively, 
and the subscript p denotes the given mixing rule. The relative deviation XiUp (further used 
as x^^^) Wcis calculated from 

rand _ ^homog 

(16) 



r)rand _ f)homc 



= : 

^hJ,P Qrand 

where QF"""^ is averaged extinction coefficient for randomly located inclusions, Q^J^°^ is 
extinction coefficient for homogeneous grains for a given mixing rule. 

The values of x[j.p cire biased by deviations related to Equation [15] in the following way: 



gqrand 



a 



rand /^homog / /^rand\—2 ^ /ITN 



An example of the dependence of x*"^"* on the wavelength for the mixture of carbon 
(matrix) and graphite (inclusions) for 20% of graphite inclusions with error bars A^*^^-* is 
shown in Figure [5l 

In order to choose the best mixing rule in the whole considered range of wavelengths 
according to the Draine criteria we have calculated the goodness-of-fit parameter x^p] (further 
used as x*-^-*) 

xS = ii:x{t (18) 

^ 1=1 

where L is the number of wavelengths for which the extinction coefficient has been calculated. 
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Figure 6. Goodness- of- fit parameter x^^^ for scattering versus volume fractions of inclusions. 
(2) 

The standard deviation of Xp j foi" ^ given jth volume fraction is calculated as: 



(1) 



^ 1 = 1 



(19) 



and is shown as error bars in Figures [9ll28l marked with letter f. 

The dependence of x^^^ on wavelength for the studied mixtures and volume fractions from 
10% to 50% with 10% step is shown in Figures marked with letters a-e. The inspection 
of the figures allows to determine the best fit of mixing rule in different wavelength ranges 
with a given volume fraction. In Figures [9ll28] marked with letter f the best fit in the whole 
range of wavelengths depending on volume fraction of inclusions has been shown. 

The similar procedure as for extinction has been carried out for the scattering and 
the scattering asymmetry parameter g =< cosO >. As an example we present the results 
of calculations of Xp'j with error bars for carbon (matrix) and graphite (inclusions) for 
scattering in Figure [6] and for asymmetry parameter in Figure [71 Compared to extinction the 
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Figure 7. Goodness-of-fit parameter x'^' for asymmetry parameter g versus volume fractions of inclusions. 

use of scattering does not improve the choice of the best mixing rule although in some, quite 
rare, cases it allows to choose better curves for some mixing rules. Furthermore, in case of the 
parameter g the obtained results are much worse. In Figure M we have shown an example 
of influence of different mixing rules on the normalised extinction E{X — V)/E{B — V) 
calculated from Mie theory. The dependence on 1/A is shown for grains with radius 0.02yum 
and consisting of ice (matrix) and 50% graphite inclusions. 

7 RESULTS AND DISCUSSION 

For the symmetrical rules studied in this paper (Lichtenecker, Bruggeman, Looyenga) the 
values of the resulting permittivity of the mixture are the same for f = 50% both when the 
material A is an inclusion in the material B and vice versa. For the values of x*-^-* and x*-^'* 
there are slight differences for f = 50%. It is caused by a different arrangement of inclusions 
when the material A is an inclusion in the material B than for an inverse situation (B - 
inclusion, A - matrix). Different arrangements of the same amount of inclusions give slightly 
different values of QIJ"^. 

7.1 Mixture of "astronomical" silicate and ice 

We have considered the " astronomical" silicate inclusions in the icy matrix. For the volume 
fraction smaller than 10% the best agreement of extinction coefficient of the mixture obtained 
from the mixing rule with the one of a mixture with "random" arrangement of inclusions in 
the whole wavelength range have been obtained for the Lichtenecker mixing rule - x^^^ has 
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Figure 8. Example of influence of the used mixing rule on normalised extinction. Ice - matrix, graphite - inclusions (50%) 



Table 1. Best mixing rules for the volume fraction of inclusions higher than 17% for the mixture of "astronomical" silicate in 
ice matrix 



Mixing rule Xmax 

Maxwell-Garnett 0.022 

Lichtenecker (symmetrical) 0.025 

Hanay 0.033 

Bruggeman (symmetrical) 0.048 

Looyenga (symmetrical) 0.058 



the smallest value. The Lichtenecker rule x*'^^ is the smallest up to 17% of volume fraction 
of inclusions. Above the 17% the best rules are listed in Table [1] and it is seen that the 
differences in x^'^^ between Lichtenecker and Maxwell-Garnett rules are very small. These 
results are displayed in Fig. O 

In case of ice inclusions in "astronomical" silicate matrix for the volume fractions of 
inclusions up to 13% the best is the Looyenga rule but not much worse is the Lichtenecker 
one - in both cases x^"^^ < 0.01. Above 13% of volume fraction of inclusions the best rules 
are listed in Table [2] and shown in Fig. [TOl 

Table 2. Best mixing rules for the volume fraction of inclusions higher than 13% for the mixture of ice in "astronomical" 
silicate matrix 



Mixing rule 


(2) 
\max 


Lichtenecker (symmetrical) 


0.022 


Looyenga (symmetrical) 


0.042 


Bruggeman (symmetrical) 


0.043 


Hanay 


0.059 


Maxwell-Garnett 


0.071 



© 2008 RAS, MNRAS 000, [T]-?? 



14 N.Maron and O. Maron 

7.2 Mixture of a - SiC and ice 

In case of a — SiC inclusions in the ice matrix (Fig. [TT!) up to 20% of volume fraction of 
inclusions the best rule is the Lichtenecker one. In the range from 20 to 40% slightly better 
from the Lichtenecker is the Maxwell- Garnett mixing rule (Ax^^) < 0.004). Above 40% agam 
the best is Lichtenecker mixing rule. 

For the inverse case (Fig. [T2|) the Looyenga rule is best for volume fraction of inclusions 
up to 12% and in the range from 12% to 50% the Lichtenecker mixing rule gives the best 
results. Next in the order of goodness are Bruggeman, Hanay and Maxwell-Garnett. The 
best mixing rule for the mixtures of a — SiC and ice is the Lichtenecker rule. For both cases 
of such mixtures the applicability of mixing rules has been examined for different wavelength 
ranges. 



7.3 Mixture of graphite and ice 

For the graphite inclusions in ice matrix up to 20% of volume fraction of inclusions the 
best results gives the Maxwell-Garnett rule and above the 20% the Hanay rule is best. 
Unfortunately both rules are asymmetrical and therefore their applicability for high values 
of volume fraction of inclusions is limited. 

In case of graphite matrix with ice inclusions in the whole volume fraction range the 
best results have been obtained for the Lichtenecker mixing rule. Figures [13] and [H] show 
the results of calculations for the above cases. 



7.4 Mixture of carbon and ice 



Up to 35% of volume fraction of carbon inclusions in the ice matrix the best is the Maxwell- 
Garnett rule and above that value the best is the Lichtenecker one. Next best rules are: 
Hanay, Bruggeman and Looyenga (Fig. [T5|) . 

In the inverse case for the whole range of volume fractions the best results gives the 



Lichtenecker rule and next best rules are: Looyenga, Hanay and Maxwel 



Draine 



■ Garne tt (Fig. [IE]). 
(Il988l ) criteria the 



For both cases of the mixture of ice and carbon considering the 
applicability of the rules has been examined for different wavelength ranges. 
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Table 3. Best mixing rules for Fe (inclusion) - ice (matrix) 



Volume fraction of inclusions [%] 


Mixing rule 


0-10 


Maxwell— Garnett 


10-33 


Hanay 


33-50 


Lichtenecker 



7.5 Mixture of Fe and ice 

Hanay (Bruggeman's asymmetric) formula corresponds best with experiment for 



ences of complex permittivities between metallic inclusion and dielectric matrix (Merill et al. 



large differ- 



Merill et al 



fll99 



19991 ). This fact is confirmed for the inclusions of iron in ice (Fig. [T711 
also stated on the basis of experiment that the Looyenga rule is only valid for low contrast 
between inclusion and matrix e^, permittivity and thus is not appropriate in the metallic 
limit which is clearly seen in Fig. [T71 The best mixing rules for this mixture are listed in 
Table El 

The mixture of ice inclusions in iron matrix satisfies the iDraind (1l988l ) criteria for the 



15% to 50% of ice content in iron matrix. In this range of volume fraction the best mixing 
rule is the formula by Lichtenecker and then down to the worst: Looyenga, Bruggeman, 
Hanay and Maxwell- Garnett (Fig. [T8|l . 



7.6 Mixture of graphite and carbon 

In the whole range of volume fractions of graphite inclusions in carbon the best rule is the 
one by Lichtenecker and then in order of decreasing goodness: Maxwell-Garnett, Hanay and 
Looyenga (Fig. [T9ll . Below 10% the quality factor for Maxwell-Garnett is slightly smaller 
than for the Lichtenecker rule. 

For carbon inclusions in graphite matrix in the whole range of volume fractions the 
best rule is the one by Lichtenecker and then in order of decreasing goodness: Looyenga, 
Bruggeman, Hanay and Maxwell-Garnett (Fig. [20]) . 

7.7 Mixture of SiC and carbon 

For the mixture of SiC inclusions in carbon matrix in the whole range of volume fractions 
of inclusions and wavelengths the best results have been obtained by using the Lichtenecker 
mixing rule. Next in decreasing order have been: Looyenga, Bruggeman, Hanay and Maxwell- 
Garnett (Fig. EH). 

For the mixture of carbon inclusions in SiC matrix the Lichtenecker mixing rule gives the 
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Table 4. The choice of best mixing rules depending on volume fraction of pores in silicate for different wavelengths 



VoluiiiG iractioii 


wavelength 


Adixing rule 


of inclusions [%] 


range [fJ.m] 






9 1 n 
u. z — lU 


Looyenga 




lU— OU 


Looyenga, Lichtenecker 




>60 


Bruggeman, Maxwell-Garnett 




9 10 


Looyenga 






Lichtenecker 




50-100 


Looyenga 




>100 


Bruggeman, Maxwell-Garnett 


25 


0.2-10 


Looyenga 




10-60 


Lichtenecker 




>60 


Looyenga 


30 


0.2-10 


Looyenga 




10-80 


Lichtenecker 




>80 


Looyenga 


40 


<10 


Looyenga 




>10 


Lichtenecker 


50 


<1.5 


Looyenga, Bruggeman 




1.5-6 


Lichtenecker 




>10 


Lichtenecker 



best results in the whole range of volume fractions of inclusions and wavelengths. Next in 
decreasing order have been: Maxwell-Garnett, Hanay, Bruggeman and Looyenga (Fig. [221) . 



7.8 Porous structures 



We have considered the " astronomical silicate" , carbon, iron and ice containing the vacuum 
pores which have been treated as inclusions of the same dimensions randomly distributed 
with refractive index m = 1.0 + il ■ 10~^°. Calculations have been carried out in the same 
way as in previous cases with 50% of inclusions (porosity). Similar calculations for porous 
materials with refractive indices characteristic for " dirty i ce" , silicate and amorp hous carbon 
in the visual wavelengths range have been carried out by 



Voshchinnikov et al. 



(120071 ). 



7.8.1 Astronomical silicate with vacuum inclusions 

The Figure [23] f shows that up to 31% of volume fraction of inclusions the best results are 
obtained from the Looyenga rule and above that value from the Lichtenecker mixing rule. 
Next in the order of goodness are Bruggeman, Hanay and Maxwell-Garnett. For different 
volume fractions of inclusions in different wavelength ranges the best mixing rules are listed 
in Table H] which is the summary of results shown in Figures [23] a-e. 
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7.8.2 SiC with vacuum inclusions 

From Figure IMf it is seen that the dependance of x*^^^ versus volume fraction of inclusions is 
similar to the same dependance for silicate with vacuum inclusions (Fig. [231F). In the range 
up to 30% the best results are obtained with Looyenga mixing rule. The other rules give 
almost the same quite good values of fit parameter x^'^^ up to 20% of volume fractions of 
inclusions. From 20% to 30% the best mixing rules are Looyenga, Lichtenecker, Bruggeman, 
Hanay, Maxwell- Garnett. Above 30% the best results are obtained from Lichtenecker rule 
and next from Looyenga, Bruggeman, Hanay and Maxwell-Garnett. The dependance of 
fitting parameter x^^^ on the wavelength for different mixing rules and volume fractions of 
inclusions are shown in Figures [2l] a-e. 

7.8.3 Graphite with vacuum inclusions 

The Fig. [25]f shows that in the range from 5% to 50% of volume fraction of inclusions (pores) 
the best mixing rule is the Lichtenecker one. Next best are Looyenga, Bruggeman, Hanay and 
Maxwell-Garnett. For different volume fractions of inclusions in various wavelength ranges 
the best mixing rules have been listed in Table [6] which summarises the results shown in 
Figures [25] a-e. 

7.8.4 Carbon with vacuum inclusions 

From the Figure [26] f one can see that in the range from 5% to 50% of volume fraction of 
inclusion (pores) the best mixing rule is the Lichtenecker one. The next best are Looyenga, 
Bruggeman, Hanay and Maxwell-Garnett. For different volume fractions of inclusions in 
different wavelength ranges the best mixing rules are listed in Table [5] which is the summary 
of figures Figures [26] a-e. 



7.8.5 Iron with vacuum inclusions 



Due to the 



Draind (119881 ) criteria the applicability of mixing rules is limited to narrow 



wavelength ranges especially for low volume fractions of inclusions. For example, the wave- 
length range for 10% of volume fraction of inclusions is for A = 0.48 ~- 0.52/xm and for 50% 
A = 0.16 ^0.85/im. 

The Figure [27] f shows that in the range from 5% to 50% of volume fractions of vacuum 
inclusions the best results have been obtained using the Lichtenecker rule. Only in case of 
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Table 5. The choice of best mixing rules depending on volume fraction of pores in carbon for different wavelengths 



VoluniG traction 


wavelengtn 


Alixing rule 


of inclusions [%] 


range [^m] 




1 n 
iU 


no n 


All rules 




0.5-11 


Looyenga 




11-150 


Lichtenecker 


on 


no ft 


Lichtenecker 




fi 10 
u — lU 


Looyenga 




10-40 


Hanay, Bruggeman, Maxwell-Garnett 




>40 


Lichtenecker 


oU 


no 1 n 


Lichtenecker 




10-150 


Looyenga 


40 


0.2-0.75 


Looyenga 




0.75-30 


Lichtenecker 




30-53 


Looyenga 




>53 


Bruggeman 


50 


<0.6 


Looyenga 




0.6-100 


Lichtenecker 




>100 


Maxwell-Garnett 



50% of volume fraction of inclusions below 0.4/im better results gives the Bruggeman rule. 
Essentially all rules give equally good results. 

7.8.6 Ice with vacuum inclusions 

From Fig. [28] f one can clearly see that in the whole range of volume fractions of inclusions 
the best results have been obtained from the Looyenga rule. Next in order of goodness are 
Bruggeman, Hanay, Lichtenecker and Maxwell-Garnett. In the range of very small values of 
k for ice (order of 10~^) only the Wiener criterion has been used because ice in this range 
acts as a very good dielectric. 

8 CONCLUSIONS 

Considering the presented results it is clearly seen that different mixing rules play important 
role depending on examined wavelengths. Nevertheless, from all chosen mixing rules for the 
considered materials in most cases (13 out of 20) the best results have been obtained using 
the Lichtenecker mixing rule. In 5 cases the Lichtenecker rule is best only for some volume 
fraction of inclusions. In case of graphite inclusions in ice the best mixing rule is the Maxwell- 
Garnett one for up to 20% of volume fraction and Hanay for higher volume fractions. The 
Looyenga mixing rule gives best results for porous ice. In case of interstellar or circumstellar 
grains the processes leading to their nucleation and growth are complicated and depending 
on not always known, changing in time and space physical and chemical conditions. It causes 
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Table 6. The choice of best mixing rules depending on volume fraction of pores in graphite for different wavelengths 



Volume fraction 


wavelength 


Mixing rule 




of inclusions [%] 


range 


pm] 






10 


0.37— 


0.44 


Lichtenecker (best) 










Looyenga, Bruggeman, Maxwell-Garnett, Hanay 






44— 




Bruggeman, Maxwell-Garnett, Hanay (almost the same) 








Looyenga, Lichtenecker. For all rules x <0.01 






0.75— 


1.85 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






1.85- 


2.30 


Looyenga (best) 










Lichtenecker, Bruggeman, Hanay, Maxwell-Garnett 




20 


n '^7 

V.O 1 — 


41 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






41— 




Bruggeman, Maxwell-Garnett, Hanay (almost the same) 








Looyenga, Lichtenecker. For all rules x <0.01 






0.55— 


2.40 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






2.40- 


2.70 


Looyenga (best), Lichtenecker 










Maxwell-Garnett, Hanay, Bruggeman (almost the same) 


30 


u.oz— 


U.OO 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






0.35- 


0.50 


Bruggeman, Maxwell-Garnett, Hanay (almost the same) 








Looyenga, Lichtenecker. For all rules x <0.02 






0.52- 


2.50 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






2.50- 


3.35 


Lichtenecker (best) 










Looyenga, Looyenga, Maxwell-Garnett, Bruggeman 




40 


0.16- 


0.45 


All rules give similar results, maximum x{^^ ~ G) = 


0.054 




0.45 


-2.5 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






2.5-4 


Lichtenecker (best) 










Looyenga, Maxwell-Garnett, Hanay, Bruggeman 




50 


0.15- 


0.36 


All rules give similar results, maximum x{^ ~ G) = 


0.081 








Lichtenecker (best) 






0.36 


-2.6 


Lichtenecker (best) 










Looyenga, Bruggeman, Hanay, Maxwell-Garnett 






2.6 


-3.7 


Lichtenecker (best) 










Looyenga, Bruggeman, Maxwell-Garnett, Hanay 






4.2 


-4.8 


Lichtenecker (best) 










Looyenga, Maxwell-Garnett, Bruggeman, Hanay 






4.8 


-5.2 


Lichtenecker (best) 





Maxwell-Garnett, Looyenga, Hanay, Bruggeman 



difficulties in determining the constituents of grains, their topology, shape and parameters 
of their size distribution. The choice of the best mixing rule decreases the number of free 
parameters when choosing the computed extinction curves to compare with the interstellar or 
circumstellar extinction. Components of the mixtures for which the best results are obtained 
from Lichtenecker rule m ay form the multicomponent grains, as for example proposed by 



Mathis &: Whiffenl (119891 ). with better derived effective dielectric function. Therefore, it is 



justified to apply the mixing rules chosen based on the numerical experiments. 
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Figure 9. Ice - matrix, silicate - inclusions (a-e volume fractions of inclusions from 10% to 50% with 10% step, f - best fit in 
the whole range of wavelengths depending on volume fraction of inclusions) © 2008 RAS, MNRAS 000, [l]-?? 
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Figure 10. Silicate - matrix, ice - inclusions (same as in Fig.|9j 
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Figure 11. Ice - matrix, SiC - inclusions (same as in Fig. |9]l 
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Figure 12. SiC - matrix, ice - inclusions (same as in Fig.[9ll 
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Figure 13. Ice - matrix, graphite - inclusions (same as in Figj9] 
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Figure 15. Carbon - matrix, ice - inclusions (same as in Fig.|9ll 

© 2008 RAS, MNRAS 000, [l]-?? 



1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 



~i — I I 1 1 1 1 1| — I — I I 1 1 1 1 1 1 — I — I I 1 1 1 1 1| 

Lichtenecker 

Bruggeman 

Maxwell-Garnett 

Looyenga | 

Hanky I \ 




J I I I I 



iL 



0.1 1 10 100 

wavelength [|a,m] 




.1 1 10 100 

wavelength [|j,m] 




1 10 100 

wavelength [|a,m] 



1.2 
1 

0.8 
0.6 
0.4 
0.2 




-0.2 



1 — I I 1 1 1 1 II — I — I I 1 1 1 1 1 1 — I — I I 1 1 1 1 1 



T 




J I I I I 




1 10 100 

wavelength [|j,m] 



0.4 



0.3 



Ci 0.2 



0.1 







1 1 1 1 1 1 1 1 1 


1 ' 1 ' 


1 1 1 1 1 1 1 


z 


\ 








"I. ~ 


X 






± ^ 










,z — 


- -p. _ 




f 


1 , 1 , 1 , 1 , 


1,1, 


,1,1,1," 



0.1 1 10 100 

wavelength [|a,m] 



5 10 15 20 25 30 35 40 45 50 55 
volume fraction of inclusions [%] 



Figure 16. Ice - matrix, carbon - inclusions (same as in Fig. |9]l 
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Figure 18. Fo - matrix, ice - inclusions (same as in Fig. |9]l 
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Figure 19. Carbon - matrix, graphite - inclusions (same as in Fig.|9ll 
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Figure 21. SiC - matrix, carbon - inclusions (same as in Fig-^^goos RAS MNRAS 000 [T]-?? 
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Figure 22. Carbon - matrix, SiC - inclusions (same as in Fig.|9]l 
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Figure 23. Silicate - matrix, vacuum - inclusions (same as in ^^-^Qg j^^g MNRAS 000 [T}^ 
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Figure 24. SiC - matrix, vacuum - inclusions (same as in Fig. |9j 
© 2008 RAS, MNRAS 000, [l]-?? 



0.1 



0.05 



0.1 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




0.1 




_i I I I I I 1 1 



J I I I I I 1 1 1 



wavelength [|a,m] 



wavelength [|j,m] 




0.5 



0.4 



0.3 



0.2 



0.1 







' I ' I ' I ' I ' I ' I ' I ' I 



I 

/I H 




5 10 15 20 25 30 35 40 45 50 55 
wavelength [|j,m] volume fraction of inclusions [%] 

Figure 25. Graphite - matrix, vacuum - inclusions (same as in ^J?J^g j^^S MNRAS 000 [l]-?? 
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